Criticality without frustration for quantum spin-1 chains 



(N 

o 

> 

O 






> 

o 

00 

rn 

O 



X 



Sergey Bravyi/ Libor Caha,^ Ramis Movassagh,'^ Daniel Nagaj,^ and Peter W. Shor^ 

^IBM Watson Research Center, Yorktown Heights NY 10598 

^Faculty of Informatics, Masaryk University, Botanick 68a, Brno, Czech Republic 

^Department of Mathematics, Massachusetts Institute of Technology, Cambridge MA 02139 

'^Research Center for Quantum Information, Slovak Academy of Sciences, Bratislava, Slovakia 

(Dated: March 1, 2013) 

Frustration-free (FF) spin chains have a property that their ground state minimizes all individual 
terms in the chain Hamiltonian. We ask how entangled the ground state of a FF quantum spin-s 
chain with nearest-neighbor interactions can be for small values of s. While FF spin-1/2 chains are 
known to have unentangled ground states, the case s = 1 remains less explored. We propose the 
first example of a FF translation- invariant spin-1 chain that has a unique highly entangled ground 
state and exhibits some signatures of a critical behavior. The ground state can be viewed as the 
uniform superposition of balanced strings of left and right parentheses separated by empty spaces. 
Entanglement entropy of one half of the chain scales as ^ logn -f 0(1), where n is the number of 
spins. We prove that the energy gap above the ground state is polynomial in 1/n. The proof relies 
on a new result concerning statistics of Dyck paths which might be of independent interest. 



The presence of long-range entanglement in the ground 
states of critical spin chains with only short-range inter- 
actions is one of the most fascinating discoveries in the 
theory of quantum phase transitions [lH3|. It can be 
quantified by the scaling law S{L) ^ logL, where S{L) 
is the entanglement entropy of a block of L spins. In 
contrast, non-critical spin chains characterized by a non- 
vanishing energy gap obey an area law [^-|6| asserting 
that S{L) has a constant upper bound independent of L. 

One can ask how stable is the long-range ground state 
entanglement against small variations of Hamiltonian pa- 
rameters? The scaling theory predicts [2, lZ| that a chain 
whose Hamiltonian is controlled by some parameter g 
follows the law S{L) ^ logi only if L does not exceed 
the correlation length ^ ^ Iff — 5c |~'', where i^ > is the 
critical exponent and gc is the critical point. For larger 
L the entropy S{L) saturates at a constant value. Hence 
achieving the scaling S{L) ^ logi requires fine-tuning of 
the parameter g with precision scaling polynomially with 
1/ L posing a serious experimental challenge. 

The stringent precision requirement described above 
can be partially avoided for spin chains described by 
frustration-free Hamiltonians. Well-known (non-critical) 
examples of such Hamiltonians are the Heisenberg fer- 
romagnetic chain [8|, the AKLT model l9|, and parent 
Hamiltonians of matrix product states [10|, |ll|. More 
generally, we consider Hamiltonians of a form H ~ 
^ ■5'jHj^j+i, where Iljj+i is a projector acting on spins 
j,j + 1 and gj > are some coefficients. The Hamil- 
tonian is called frustration-free (FF) if the projectors 
Hj j+i have a common zero eigenvector -0. Such zero 
eigenvectors tp span the ground subspace of H. Clearly, 
the ground subspace does not depend on the coefficients 
gj as long as they remain positive. This inherent sta- 
bility against variations of the Hamiltonian parameters 
motivates a question of whether FF Hamiltonians can 
describe critical spin chains. 



In this Letter we propose a toy model describing a 
FF translation- invariant spin-1 chain with open bound- 
ary conditions that has a unique ground state with a loga- 
rithmic scaling of entanglement entropy and a polynomial 
energy gap. Thus our FF model reproduces some of the 
main signatures of critical spin chains. In contrast, it was 
recently shown by Chen et al [1^ that any FF spin-1/2 
chain has an unentangled ground state. Our work may 
also offer valuable insights for the problem of realizing 
long-range entanglement in open quantum systems with 
an engineered dissipation. Indeed, it was shown by Kraus 
et al [13[ and Verstraete et al IJ] that the ground state of 
a FF Hamiltonian can be represented as a unique steady 
state of a dissipative process described by the Lindblad 
equation with local quantum jump operators. A proposal 
for realizing such dissipative processes in cold atom sys- 
tems has been made by Diehl et al 15 1. 

Main results. We begin by describing the ground 
state of our model. The three basis states of a single spin 
will be identified with a left bracket I = [, right bracket 
r =], and an empty space represented by 0. Hence a 
state of a single spin can be written as a|0) -I- /3|Z) -(- 7|r) 
for some complex coefficients a, /3, 7. For a chain of n 
spins, basis states \s) correspond to strings s S {0, 1, r}". 
A string s is called a Motzkin path [la] iff (i) any initial 
segment of s contains at least as many I's as r's, and (ii) 
the total number of I's is equal to the total number of r's. 
For example, a string lllrOrlOrr is a Motzkin path while 
lOlrrrllr is not since its initial segment lOlrrr has more 
r's than Vs. By ignoring all O's one can view Motzkin 
paths as balanced strings of left and right brackets. We 
shall be interested in the Motzkin state \A4n) which is the 
uniform superposition of all Motzkin paths of length n. 
For example, jA^z) -- |00) + \lr), {Ma) - |000) + \lrO) + 
I /Or) + \Olr), and 



IMa) ^ 10000) + \00lr) + \OlOr) + \lOOr) 
+ \OlrO) + \lOrO) + |/rOO) + |//rr) 



+ \lrlr 



Let us first ask how entangled is the Motzkin state. For 
a contiguous block of spins A, let pA = Tr^^^ |A^n)(A^n| 
be the reduced density matrix of A. Two important mea- 
sures of entanglement are the Schmidt rank x{A) equal 
to the number of non-zero eigenvalues of pA, and the 
entanglement entropy S{A) = — Tr p^ logj pa ■ We will 
choose A as the left half of the chain, A = {1, . . . , n/2}. 
We show that 

X{A) ^l + n/2 and S{A) = - log2 n + c„ (1) 

where lim„_i.oo c„ — 0.14(5). The linear scaling of the 
Schmidt rank stems from the presence of locally un- 
matched left brackets in A whose matching right brackets 
belong to the complementary region B — [l,n\\A. The 
number of the locally unmatched brackets m can vary 
from to n/2 and must be the same in A and B lead- 
ing to long-range entanglement between the two halves 
of the chain. 

Although the definition of Motzkin paths may seem 
very non-local, we will show that the state \Mn) can 
be specified by imposing local constraints on nearest- 
neighbor spins. Let H be a projector onto the three- 
dimensional subspace of C^ ^ C^ spanned by states 
\0l) -\10), I Or) - |rO), and |00) -\lr). Our main result is 
the following. 

Theorem 1. The Motzkin state |A^„) is a unique ground 
state with zero energy of a frustration- free Hamiltonian 



log A 



iJ=|r)(r|i + |0(/|„ + 5]n,.,+i, 



(2) 



where subscripts indicate spins acted upon by a projector. 
The spectral aavil 21 of H scales polynomially with 1/n. 

The theorem remains true if H is modified by introduc- 
ing arbitrary weights gj > 1 for every projector in Eq. ([2|) . 
A polynomial lower bound on the spectral gap of H is, 
by far, the most difficult part of Theorem [1] Our proof 
consists of several steps. First, we use a perturbation 
theory to relate the spectrum of H to the one of an effec- 
tive Hamiltonian iJeff acting on Dyck paths — balanced 



strings of left and right brackets [18| . This step involves 
successive applications of the Projection Lemma due to 
Kempe et al [19|. Secondly, we map H^s to a stochas- 
tic matrix P describing a random walk on Dyck paths 
in which transitions correspond to insertions/removals of 
consecutive Ir pairs. The key step of the proof is to show 
that the random walk on Dyck paths is rapidly mixing. 
Our method of proving the desired rapid mixing prop- 
erty employs the polyhedral description of matchings in 
bipartite graphs [20[. This method appears to be new 
and might be interesting on its own right. Exact diag- 
onalization performed for short chains suggests that the 
spectral gap of H scales as A ^ 1/"'^, see Fig. [1] Our 
proof gives an upper bound A — 0{n^^''^) and a lower 
bound A = il{n~'') for some c^ 1. 



log n 



FIG. 1. The spectral gap A of the Hamiltonian H defined in 
Eq. ([21 for 3 < n < 13 obtained by the exact diagonalization. 
The dashed line shows a linear fit log A = 0.68 — 2.91 log n. 
Our numerics suggests that the first excited state of H be- 
longs to the subspace spanned by strings with exactly one 
unmatched bracket. 



Previous work. Examples of spin chain Hamiltoni- 
ans with highly entangled ground states have been con- 
structed by Gottesman and Hastings 
for local dimension d = 9 and d 
(here and below d = 2s + 1). These models exhibit 
a linear scaling of the entropy S{L) for some blocks 
of spins while the spectral gap is polynomial in 1/n. 



21|, and Irani [221 
= 21 respectively 



The model found in [2^ is FF and translation-invariant. 
Ref. [23j focused on 'generic' spin chains with a Hamil- 
tonian H = ^ ■ IljJ-^-l where the projectors Hjj+i are 
chosen randomly with a fixed rank r ^24J. The authors 
of [23| identified three important regimes: (i) frustrated 
chains, r > (i^/4, (ii) FF chains, d < r < d^/4, and (in) 
FF chains with product ground states, r < d. It was 
conjectured in [23[ that generic FF chains in the regime 
d < r < d^/A have only highly entangled ground states 
with probability one. This regime however requires lo- 
cal dimension d > 4. The new model based on Motzkin 
paths corresponds to the case d = r = 3 (ignoring the 
boundary terms) and thus it can be frustrated by arbi- 
trarily small deformations of the projectors making them 
generic. In addition, results of ,23.] imply that examples 
of FF spin-1 chains with highly entangled ground states 
have measure zero in the parameter space. The question 
of whether matrix product states specified by FF par- 
ent Hamiltonians can exhibit quantum phase transitions 
has been studied by Wolf et al [25]. However, the mod- 
els studied in J25| have bounded entanglement entropy, 
5(L) = 0(1). 

Hamiltonian. Let us now construct a FF Hamilto- 
nian H whose unique ground state is \Ain)- First we need 
to find a more local description of Motzkin paths. Let 
S = {0, 1, r}. We will say that a pair of strings s,t ^ 'S"' 
is equivalent, s ~ i, if s can be obtained from t by a 



sequence of local moves 

00 ^ Ir, 01 <-> 10, Or ^ rO. 



(3) 



These moves can be applied to any consecutive pair of 
letters. For any integers p,q > such that p + q < n 
define a string 



^p,q 



r...rO...Ol...l. 

p n—p—q q 



Lemma 1. Any string s £ E" is equivalent to one and 
only one string Cp^q. A string s £ S" is a Motzkin path 
iff it is equivalent to the all-zeros string, s ^ co,o- 

Proof. Indeed, applying the local moves Eq. (jH]) one 
can make sure that s does not contain substrings Ir or 
10 . . . Or. If this is the case and s contains at least one I, 
then all letters to the right of I are I or 0. Similarly, if 
s contains at least one r, then all letters to the left of r 
are r or 0. Since we can swap with any other letter by 
the local moves, s is equivalent to Cp^q for some p,q. It 
remains to show that different strings Cp^q are not equiv- 
alent to each other. Let Lj{s) and Rj{s) be the number 
of I's and r's among the first j letters of s. Suppose 
Cp^q ~ Cp'^q' and p > p'. Then Rp{s) — Lp{s) < p' for 
any string s equivalent to Cpi^qi. This is a contradiction 
unless p — p' . Similarly one shows that q — q' . D 

The lemma shows that the set of all strings E" can be 
partitioned into equivalence classes Cp^q, such that Cp^q 
includes all strings equivalent to Cp_q. In other words, 
s € Cp,q iff s has p unmacthed right brackets and q un- 
matched left brackets. Accordingly, the set of Motzkin 
paths jMn coincides with the equivalence class Co,o- 

Let us now define projectors 'implementing' the local 
moves in Eq. ([3]). Define normalized states 

|0)^|OO)-|/r), |^V|0/)-|/0), \r)-\Or}-\rO) 

and a projector H = \(f)){<f)\ + |?A0(V''I + IV''')(V'1- Ap- 
plication of n to a pair of spins j, j -I- 1 will be denoted 
IIjj+i. If some state "0 is annihilated by every projector 
IIjj+i, it must have the same amplitude on any pair 
of equivalent strings, that is, {s\ijj) = (t|'0) whenever 



t. It follows that a Hamiltonian iJ^ 



E"",'n 



j=i iij,i-(-i 



is FF and the ground subspace of 7J^ is spanned by pair- 
wise orthogonal states \Cp^q), where \Cp^q) is the uni- 
form superposition of all strings in Cp^q. The desired 
Motzkin state \Mn) — |C'o,o) is thus a ground state of 
H^. (It is worth mentioning that not all states \Cp^q) 
are highly entangled. For example, |C„_o) = |r)®" is 
a product state.) How can we exclude the unwanted 
ground states \Cp^q) with p ^ oi q ^ 07 We note 
that Co.o is the only class in which strings never start 
from r and never end with I. Hence a modified Hamil- 
tonian H = |?')(r|i -I- |0(^|n + Hr^ that penalizes strings 
starting from r or ending with / has a unique ground 



state |Co,o)- This proves the first part of Theorem [T] 
We can also consider weighted Hamiltonians Hr^{g) = 
Y.]Zl 9j^3,j+i and H{g) = go\r){r\i + gn\l){l\n + H^{g), 
where go, ■ • ■ ,5n > 1 are arbitrary coefficients. One can 
easily check that the ground state of H{g) does not de- 
pend on g and H{g) > H . It implies that the spectral 
gap of H{g) is lower bounded by the one of H. 

Entanglement entropy. We can now construct the 
Schmidt decomposition of the Motzkin state. Let A = 
{1, . . . , n/2} and B = {n/2 + I, . . . , n} be the two halves 
of the chain (we assume that n is even) . For any string 
s e E" let SA and sb be the restrictions of s onto A 
and B. We claim that s is a Motzkin path iff sa ~ co,m 
and sb ~ Cm,,o for some < to < n/2. Indeed, sa 
[sb) cannot have unmatched right (left) brackets, while 
each unmatched left bracket in sa must be matched with 
some unmatched right bracket in sb- It follows that the 
Schmidt decomposition of \M.n) can be written as 



n/2 
\Mn) = E 



/Pm IC'o.m)^ ® |C'm,o)B, 



(4) 



where \Cp^q) is the normalized uniform superposition of 
all strings in Cp,q and Pm are the Schmidt coefficients 
defined by 

Pti 



|g0,rn(n/2)|- 

|C^o,o(n)| 



(5) 



Here we added an explicit dependence of the classes Cp^q 
on n. For large n and m one can use an approximation 
Pm ^ rn? exp {—Zm^/n), see the Supplementary Material 
for the proof. Note that Pm achieves its maximum at 
TO* « \/n/i. Approximating the sum X^mP™ = 1 by 



an integral over a = m/^Jn one gets p„ 
where qa is a normalized pdf defined as 



-1/2 



qa(r. 



qa = Z-'a^e--"^ , Z= daa^e 



,2„-3a^ 



4 •33/2' 



It gives 

/•oo 

S'(A) := - ^ Pm log2 Pm ~logy/n- da q^ log2 qa ■ 

Evaluating the integral over a yields Eq. (jlj. The 
approximation pm ~ 'n-"^ qa{m) also implies that 
maxmPm = 0{n^^''^). This bound will be used below 
in our spectral gap analysis. We conjecture that one can 
achieve a power law scaling of S{A) in Eq. ([T]) by intro- 
ducing two types of brackets, say I = [, r =], V = {, 
and r' =}, such that bracket pairs Ir and I'r' are cre- 
ated from the 'vacuum' 00 in a maximally entangled state 
{\lr) + \l'r'))/^/2. The local moves Eq. ^ must be mod- 
ified as Ox -ir^ xO, where x can be either of l,r,l' ,r' , and 
00 <-^ {Ir + l'r')/\/2. We expect the modified model with 
two types of brackets to obey a scaling S{A) ~ y/n, while 



its gap will remain lower bounded by an inverse polyno- 
mial. 

Spectral gap: upper bound. Let A2 > be the 

smallest non-zero eigenvalue of the Hamiltonian defined 
in Eq. ^. We shall use the fact that the ground state 
\^An) is highly entangled to prove an upper bound A2 < 
0(n~^/^). Fix any k £ [0,n/2] and define a 'twisted' 
version of the ground state: 



n/2 
m=0 



/Prn\Co,m)A €5 |Cm,o)i 



where 9„i = for < tti < fc and 0m = tt otherwise. 
Note that \(p) and \Mn) have the same reduced den- 
sity matrices on A and B. Hence {(l)\H\ip) — (0|nctif|0), 
where Ucut = n„/2,n/2+i- Since laaxmPm = 0(n"i/2) 
and J2mP"i ~ 1' there must exist k £ [0, n/2] such that 
So<m<fcP'" = f /2 ± e for some e = 0{n-^^^). This 
choice of k ensures that (0|A^„) ~ '^mP^nS-^^"^ < 2e, 
that is is almost orthogonal to the ground state. De- 
fine a normalized state |0) ^ If/*) ~ (■^ri|0) • \M.n)- Then 



n. 



< 



n. 



((/.|7W„) = and (0|i?|0) = ,,, _ ,, , 

0(e). The difference ((/)|nc„t|(/)) - {Mn\^cut\Mn) gets 
contributions only from the terms ?Ti = fc, fc ± 1 in the 
Schmidt decomposition, since Ilcut can change the num- 
ber of unmatched brackets in A and B at most by one. 
Since {Mn\^cut\Mn) = 0, we get 



mr 



< 0{pk+Pk-i +Pk+i) = 0{n-^/^). 



0{n 



-1/2N 



Therefore A2 is at 



We arrive at {(f>\H\(f> 
most 0(n~^/^). 

Spectral gap: lower hound. It remains to prove a 
lower bound A2 > n^'~'^^\ Let Tip.q be the subspace 
spanned by strings s G Cp^q and Hm = 'Ho,o be the 
Motzkin space spanned by Motzkin paths. Note that 
H preserves any subspace 'Hp,q and the unique ground 
state of H belongs to V-m ■ Therefore it suffices to derive 
a lower bound n"^^^^ for two quantities: (i) the gap of H 
inside the Motzkin space %m^ and (ii) the ground state 
energy of H inside any 'unbalanced' subspace T-Lp^q with 
p 7^ or q 7^ 0. Below we shall focus on part (i) since 
it allows us to introduce all essential ideas. The proof of 
part (ii) can be found in the Supplementary Material. 

Recall that a string s G {l,r}'^"^ is called a Dyck path 
iff any initial segment of s contains at least as many Vs 
as r's, and the total number of Vs is equal to the to- 
tal number of r's. For example, Dyck paths of length 6 
are lllrrr, llrlrr., llrrlr, Irlrlr, and Irllrr. The proof of 
part (i) consists of the following steps: 
Step 1. Map the original Hamiltonian H acting on 
Motzkin paths to an effective Hamiltonian H^s acting 
on Dyck paths using perturbation theory. 
Step 2. Map Hgg to a stochastic matrix P describing a 
random walk on Dyck paths in which transitions corre- 
spond to insertions or removals of consecutive Ir pairs. 



Step 3. Bound the spectral gap of P using the canonical 
paths method [26| . 

To construct a good family of canonical paths in Step 3 
we will organize Dyck paths into a rooted tree in which 
level-TO nodes represent Dyck paths of length 2m, edges 
correspond to insertion of Ir pairs, and each node has 
at most four children. Existence of such a tree will be 
proved using the fractional matching method [20|. 

Let I?m be the set of Dyck paths of length 2m, V 
be the union of all "Dm with 2m < n, and A^„ be the 
set of Motzkin paths of length n. Define a Dyck space 
Hd whose basis vectors are Dyck paths s £ T>. Given a 
Motzkin path u with 2m brackets, let Dyck(M) £ P^ be 
the Dyck path obtained from u by removing zeros. We 
shall use an embedding V : Hd ^ Hm defined as 



V\s)^ 



1 



y V2m 



E 

UGMn 

Dyck(u) — s 



s€Vnv„ 



One can easily check that V'^V = I, that is, V is an isom- 
etry For any Hamiltonian H, let \2{H) be the second 
smallest eigenvalue of H . 

Step 1. The restriction of the Hamiltonian Eq. ([2]) onto 
the Motzkin space T-Lm can be written as H — Hmove + 
Hint , where Hmove describes freely moving left and right 
brackets, while Hint is an 'interaction term' responsible 
for pairs creation. More formally, H, 



and Hi, 



'-move 
jmove 



En—l Y\move 
.7 = 1 ^^" 



^J = i " Jj' + l 

Hnt - T,'j=i n}"^!, where n'"°"'= projects onto the 
subspace spanned by \0l) — \IQ) and |0r) — |rO), while 
n"'* projects onto the state |00) — \lr). Note that the 
boundary terms in Eq. ([2]) vanish on Hm ■ We shall treat 
Hint as a small perturbation of Hmove- To this end de- 
fine a modified FF Hamiltonian H^ = Hmove + eHint, 
where < e < 1 will be chosen later. One can easily 
check that \Mn) is the unique ground state of H^ and 
^2{H) > X2{He) (use the operator inequality H > H^). 
Note that Hmove4' = iff "0 is symmetric under the moves 
01 <->■ 10 and Or <->■ rO. It follows that the ground sub- 
space of Hmove is Spanned by states V \s) with s £ V. 
To compute the spectrum of Hmove i we can ignore the 
difference between Vs and r's since Hmove is only capable 
of swapping zeros with non-zero letters. It follows that 
the spectrum of Hmove must coincide with the spectrum 
of the Heisenberg ferromagnetic spin-I/2 chain, that is, 
Y^move ^gj^ ^g replaced by the projector onto the singlet 
state |01) — 1 10), where 1 represents either I or r. Using 
the exact formula for the spectral gap of the Heisenberg 
chain found by Koma and Nachtergaele [8| we arrive at 

X2{Hmove) = 1 - COS (^) = ^{71-"^). Let 
TJeff = V^H.ntV 

be the first-order effective Hamiltonian acting on the 
Dyck space Hd- Applying the Projection Lemma of [19[ 
to the orthogonal complement of |A^„) in T-Lm we infer 



that 



A2(i/e)>eA2(i/cff) 



Oie^)m, 



^2\Hmove) ~ 2e||iJ.j, 



(6) 



Choosing e <C n '^ guarantees that e||iJmt|| is small com- 
pared with \2{Hmove)- For this choice of e one gets 



2„4\ 



\2{H) > \2{H,) > eA2(i/off) - Oie'ri 



(7) 



Hence it suffices to prove that A2(i?off) > n^'~'^^\ 
Step 2. Recall that Hctf — V'^HintV acts on the Dyck 
space "Hd- Its unique ground state \'D) G %£> can be 
found by solving \Mn) — V\'D). It yields 



\V) 



V\M^\ 




(8) 



'2in<n 



Let 7r(s) = (s|2?)^ be the induced probability distribution 
on v. Given a pair of Dyck paths s,t G V, define 



n W n[s) 



(9) 



We claim that P describes a random walk on the set 
of Dyck paths T) such that P{s,t) is a transition prob- 
ability from s to t, and 7r(s) is the unique steady state 
of P. Indeed, since \/n{t) is a zero eigenvector of iJcff, 
one has Et^(s,0 = 1 and E^ 7^(s)-P(s>i) = 7r(i). Off- 
diagonal matrix elements {s\Hcff\t) get contributions only 
from terms — |00)(^r| and — |^r)(00| in Hint- It implies 
that {s\Hcs\t) < for s 7^ t and hence P{s,t) > 0. Fur- 
thermore, P{s,s) > 1/2 since (s|_ffcff|s) < n/2. In the 
Supplementary Material we shall prove the following. 

Lemma 2. Let s,t G V be any Dyck paths such that t 
can be obtained from s by adding or removing a single Ir 
pair. Then P(s,t) = fl{l/n^). Otherwise P{s,t) = 0. 

Let A2(P) be the second largest eigenvalue of P. From 
Eq. ^ one gets A2(i?off) = n{l - A2(P)). Hence it suf- 
fices to prove that the random walk P has a polynomial 
spectral gap, that is, 1 — \2{P) > n~'^(^\ 
Step 3. Lemma [2] tells us that P describes a random 
walk on a graph G = (2?, E) where two Dyck paths are 
connected by an edge, {s^t) G -E, iff s and t are related 
by insertion/removal of a single Ir pair. To bound the 
spectral gap of P we shall connect any pair of Dyck paths 
s,t G I? by a canonical path j{s,t), that is, a sequence 
So, si, . . . , s/ G I? such that sq — s, si ~ t, and (si, s^+i) G 
E for all i. The canonical paths theorem [2(| shows that 
1 — A2(P) > l/{pl), where I is the maximum length of a 
canonical path and p is the maximum edge load defined 



p = max , . „, — "Y 7 

{a,b)£E Tr{a)P{a,b) , , , , 



7r(s)7r(i). (10) 



The key new result that allows us to choose a good family 
of canonical paths is the following. 



Lemma 3. Let Vk be the set of Dyck paths of length 2k. 
For any k > 1 there exists a map f : T>k -^ "Dk-i such 
that (i) the image of any path s G Vk can be obtained from 
s by removing a single Ir pair, (ii) any path t G 2?fe-i has 
at least one pre-image in T^k, and (Hi) any path t G I^fc-i 
has at most four pre-images in Vk. 

The lemma allows one to organize the set of all Dyck 
paths V into a supertree T such that the root of T rep- 
resents the empty path and such that children of any 
node s are elements of f^^{s). The properties of / im- 
ply that Dyck paths of length 2m coincide with level-m 
nodes of T, any step away from the root on T corre- 
sponds to insertion of a single Ir pair, and any node of T 
has at most four children. Hence the lemma provides a 
recipe for growing long Dyck paths from short ones with- 
out overusing any intermediate Dyck paths. It should be 
noted that restricting the maximum number of children 
to four is optimal since |I?fe| — Ck ~ A^ / ^/TTk^/"^, where 
Cfc is the /c-th Catalan number. Our proof of Lemma [3] 
based on the fractional matching method can be found 
in the Supplementary Material. Five lowest levels of the 
supertree T are shown on Fig. [51 



We can now define the canonical path 7(5, i) from 
s G Vm to i G Pfc. Any intermediate state in 7(5, t) 
will be represented as uv where u G T^v is an ancestor of 
s in the supertree and v G Di" is an ancestor of t. The 
canonical path starts from u = s, u = and alternates 
between shrinking u and growing v by making steps to- 
wards the root (shrink) and away from the root (grow) 
on the supertree. The path terminates as soon as m = 
and V = t. The shrinking steps are skipped whenever 
u = 0, while the growing steps are skipped whenever 
V — t. Note that any intermediate state uv obeys 



min(|s|, \t\) < \u\ + \v\ < max(|s|, \t\). 



fll) 



Since any path 7(5, t) has length at most 2n, it suffices to 
bound the maximum edge load p. Fix the edge (a, b) G E 
with the maximum load. Let p{m, k, I', I") be the contri- 
bution to p that comes from canonical paths 7(5, i) such 
that a — uv (^ Vii^in, where 



seVr, 



teVk, ueVi,, veVi. 



and such that b is obtained from a by growing v (the case 
when b is obtained from a by shrinking u is analogous). 
The number of possible source strings s G I?m contribut- 
ing to p{m, k, I', I") is at most 4™~' since s must be a 
descendant of u on the supertree. The number of possi- 
ble target strings t G Vk contributing to p{m, k, I', I") is 
at most 4^^"' since t must be a descendant of v on the 
supertree. Taking into account that 7r(s) and Tr{t) are 
the same for all s G Vm and t G Vk we arrive at 



p(m,fc,r,/") <4" 



-k-i'-i" 7r(s)7r(t) 



TTrnTTfe 



Ti{a)P{a,b) ■Ki,+i»P{a,b) 
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FIG. 2. (Color Online) Five lowest levels of the supertree T. Nodes of T are Dyck paths — balanced strings of left and right 
brackets. Depth-fe nodes are in one-to-one correspondence with Dyck paths of length 2fc (the set Ofe). Any step towards the 
root requires removal of a consecutive [] pair. Any node has at most four children. The supertree can be described a family 
of maps / : Vt -^ "Dk-i such that /(s) is the parent of s. The maps / are defined inductively such that /([X]) = [f{X)], 
f{[\Y) = []f{Y) for any node, /([X]V) = [f{X)]Y for black nodes, and f{[X]Y) = [X]f{Y) for red (shaded) nodes. See the 
proof of Lemma |3] in the Supplementary Material for more details. 



with TTi — 4'(2;)/|A^„|. Here we used the identity 
Ti{w) — (w|2?)^ and Eq. ([8]). Lemma [2] implies that 
1/P(a, h) < n^'-^\ Furthermore, the fraction of Motzkin 
paths of length n that have exactly 21 brackets is ct; = 
Ci{^^)/\Mn\- However Q « A^/^l^/^ coincides with 4' 
modulo factors polynomial in 1/n. Hence 



p{m, k,l ,1 ) < n 



0(1) O'rnO'k 



(^I'+l" 



By definition, ct/ < 1 for all I. Also, one can easily 
check that ai as a function of I has a unique maximum at 
I » n/3 and decays monotonically away from the max- 
imum. Consider two cases. Case (1): V + I" is on the 
left from the maximum of cr;. From Eq. ([TT|) one gets 
min(TO, fc) < V + 1" and thus am(:^k < o-mm(m,fc) < f^i'+i"- 
Case (2): I' + I" is on the right from the maximum of 
CT;. From Eq. (fTT|) one gets max(m, fc) > I' + I" and 
thus amCTk < cr^^^(^rn,k) < c^i'+i"- In both cases we get a 
bound p{m, k, V , I") < n^^^'> . Since the number of choices 
for m, k, I', I" is at most n'-"-^^ , we conclude that p < n'^'^^^ 
and thus 1 - A2(P) > n-c(i). 

Open problems. Our work raises several questions. 
First, one can ask what is the upper bound on the ground 
state entanglement of FF spin-1 chains and whether the 
Motzkin state achieves this bound. For example, if the 
Schmidt rank xi^) for ^ block of L spins can only grow 
polynomially with L, as it is the case for the Motzkin 
state, ground states of FF spin-1 chains could be effi- 
ciently represented by Matrix Product States [27| (al- 
though finding such representation might be a computa- 
tionally hard problem ^2^]). One drawback of the model 
based on Motzkin paths is the need for boundary condi- 
tions and the lack of the thermodynamic limit. It would 
be interesting to find examples of FF spin- 1 chains with 
highly entangled ground states that are free from this 
drawback. We also leave open the question of whether 
our model can indeed be regarded as a critical spin chain 



in the sense that its continuous limit can be described by 
a conformal field theory. Finally, an intriguing open ques- 
tion is whether long-range ground state entanglement (or 
steady-state entanglement in the case of dissipative pro- 
cesses) in ID spin chains can be stable against more gen- 
eral local perturbations, such as external magnetic fields. 



SUPPLEMENTARY MATERIAL 
Schmidt coefficients of the Motzkin state 

In this section we compute the Schmidt coefficients Pm 
defined in Eq. ([5]) and show that for large n and m one 
can use an approximation 



Pm ^ Tn exp(— 3m- /n). 



(12) 



Let Vn^k C {Z,r}^"+'^ be the set of balanced strings of 
left and right brackets of length 2n -I- fc with fc extra left 
brackets. More formally, s G Vn^k ifi' any initial segment 
of s contains at least as many Vs as r's, and the total 
number of Vs is equal to n + k. 



Lemma 4 (Andre's reflection method). 

number of strings in 'Dn,k is 



The total 



Dn,k — 



fc + 1 



2n + fc 
1 I n 



Proof. For any bracket string s let L{s) and R{s) be the 
number of left and right brackets in s. Any s £ {I, r}2"+'= 
such that s ^ T^n^k can be uniquely represented as s = 
urv, where r corresponds to the first unmatched right 
bracket in s, while u is a balanced string (Dyck path). 
Let v' be a string obtained from w by a reflection r -H- I 
and s' = urv' . Then 



L{s') = L{u) + Liy') = R{u) + R{v) = R{s) - 1 



1 



and 



Proof of Lemma [2] 



R{s') = R{u) + 1 + R{v') = L{u) + 1 + L{v) = L{s) + 1 
= n + fc + 1. 

Furthermore, any string s' with n — \ left brackets and 
n -\- k -\- 1 right bracket can be uniquely represented as 
s' = urv' . Hence the number of strings in Vn.k is 

2n + k\ /2n + fc\ k + l /2n + k 



n — 1 J n + k + 1 \ n 



U 



One can easily check that \Cp^q{n)\ = |Co,p+g(ri)| since 
the identity of unmatched brackets does not matter for 
the counting. Let 

Mn,rn = |Co,m(n)|. 

Lemma m implies that 



Mn,^ = Y^ 



m + 1 / n \ (2i + m 



i>0 
2i-\-m<n 



l + 7TT,+ lV2z + ?71 



It can be rewritten as 



(13) 



i>0 
2i+m<n 



where 



M„ 



(m + 1) • n! 



(i + m + l)!i!(n — 2i — m)\ 



(14) 



Let a = mj \pn and f3 — {i — n/3)/y/n. Using Stirling's 
formula one can get 



M„, 



^^ 3"+iaexp(-3a'-9a/3-9/32). (15) 



27rn3/ 



We approximate the sum in Eq. P^ by integrating over 
i. Since i = ^ + [iy/n, we get di — ^/ndf3, Since the 
maximum is near i = ^, we can turn this sum into an 
integral from — oo to oo. The integral we need to evaluate 
is thus 



M„ 



-^^3"+ia / e-^^^-s^/'-^^^y^d/? 



27rn3/2 



^^3"+ia / e-9(^-"/2)^-i"^V^d/3 



27rn 



2y/Trn 
Recalling that 



^^ 3"+^ae-t"' -mexp(-3mV4n). 



_ |Co,^(n/2)p 



(16) 



we arrive at Eq. ([T 



Suppose s ^Vk and i £ 2?fe±i. Using the definition of 



P{s,t) one can easily get 



Pis,t) 



1 /ti 



77 V2A: 



E 



(M|iJ.„t|w) 



ueM„\s] veMnlt] 



Here A^„[s] = {u G Mn ■ Dyck(u) = s}. Note that 
{u\Hint\v) = —1/2 if u and w differ exactly at two con- 
secutive positions where u and v contain 00 and Ir re- 
spectively or vice versa. In all other cases one has 
{u\H„,t\v) ^ 0. 

Suppose t e I'fc+i and P{s, t) > 0. Let us fix some j g 
[0, 2k] such that t can be obtained from s by inserting a 
pair Ir between Sj and Sj+i. For any string u G A^„[s] in 
which Sj and Sj+i are separated by at least two zeros one 
can find at least one v € Ain[t] such that {u\Hint\v) = 
— 1/2. The fraction of strings u £ A4n[s] in which Sj and 
Sj+i are separated by two or more zeros is at least 1/n^ 
which implies 



Pis,t)> 



2^3 ■ 



Suppose now that t G 2?fc-i • Let us fix some j G [1, 2fc— 1] 
such that t can be obtained from s by removing the pair 
SjSj+i = Ir. For any string u G jV{„[s] in which Sj and 
Sj+i are not separated by zeros one can find at least one 
V G A^„[i] such that {u\Hint\v) = —1/2. The fraction of 
strings u G A^„[s] in which Sj and s^+i are not separated 
by zeros is at least \/n which implies 



^(^-i)>. .■ 



1 
2~2 



Proof of Lemma [3] 

Let us first prove a simple result concerning fractional 
matchings. Consider a bipartite graph G = {AU B,E). 
Let X = {xe}eeE be a vector of real variables associated 
with edges of the graph. For any vertex u let S{u) be the 
set of edges incident to u. Define a matching polytope (20| 

V = {x : Xe > for all e e E, 

1 < X! ^^ - 4' E ^« = 1 

eei5(a) eeS(b) 

for aU a G A and b G B}. 

Lemma 5. Suppose V is non-empty. Then there exists a 
map f : B ^ A such that (i) f{b) = a implies (a, b) G E, 
(ii) any vertex a £ A has at least one pre-image in B, 
and (Hi) any vertex a £ A has at most four pre-images 
in B. 

Proof. Since V is non-empty, it must have at least one 
extremal point x* G V. Let £'* C £' be the set of edges 



such that a;* > 0. We claim that E* is a forest (disjoint 
union of trees). Indeed, suppose E* contains a cycle Z 
(a closed path). Then a;* j, < 1 for all (a, 6) £ Z since 
otherwise the cycle would terminate at b. Hence < 
x* < 1 for all e G Z. Since the graph is bipartite, one 
can label edges of Z as even and odd in alternating order. 
There exists e ^ such that x* can be shifted by ±e on 
even and odd edges of Z respectively without leaving V . 
Hence x* is a convex mixture of two distinct vectors from 
V. This is a contradiction since x* is an extreme point. 
Hence E* contains no cycles, that is, E* is a forest. We 
claim that x* = 1 for all e e E* . Indeed, let T C £:* be 
the subset of edges with < a;* < 1. Obviously, T itself 
is a forest. Degree- 1 nodes of T must be in A and there 
must exist a path 7 C T starting and ending at degree-1 
nodes u,u' E A. Since < x* < 1 for all e € 7, there 
exists e 7^ such that x* can be shifted by ±e on even 
and odd edges of 7 respectively without leaving V. This 
is a contradiction since x* is an extreme point. Hence 
x* = 1 for all e e E*. We conclude that a;* G {0, 1} for 
all edges of G. The desired map can now be defined as 



fib) ^aiSx: 



1. 



D 



We are interested in the case where 

A = Vn-i and B = V^ 

are Dyck paths of semilength n — 1 and n respectively. 
Paths a e 'Dn-i and b e Vn are connected by an edge iff 
a can be obtained from b by removing a single Ir pair. 
Our goal is to construct a map / : I?„ — > Vn-i with 
the properties (i),(ii),(iii) stated in Lemma [51 Accord- 
ing to the lemma, it suffices to choose / as a stochastic 
map. Namely, for any b S I?„ we shall define a random 
variable /(&) e I'n-i with some normalized probability 
distribution. It suffices to satisfy two conditions: 

Pr[/(6) = a] > only if a can be obtained from b 
by removing a single Ir pair, 

(17) 
and 



bev. 



Pr[/(6) =a]= X„ for all a e V„ 



(18) 



Here 1 < Xn < 4 is some function of n that we shall 
choose later. We shall define / using induction in n. 

Lemma 6. For every n > 1 there exists 1 < Xn < 4 and 
a stochastic map f : !)„ — > I'ri-i satisfying Eqs. I^17]18\} . 

Proof. Any Dyck path b G 2?„ can be uniquely repre- 
sented as b = Isrt for some s € Vi, t € Vn-i-i, and 
J G [0, n — 1]. We shall specify the map / : 2?„ -^ 2?n-i 
by the following rules: 



beVn 




fib) e v„-i 


probability 


Isrt, s EVi, 1 < i < n - 


-2 


lf{s)rt 


Pi 


Isrt, s EVi, 1 < i < n- 


-2 


Isrfit) 


l^P^ 


Irt, t e 2?„_i 


Irfit) 


1 


Isr, S e Pn-l 


lf{s)r 


1 



Here we assumed that / has been already defined for 
strings of semilength up to n — 1 such that Eqs. (|17ll8p 
are satisfied. By abuse of notation, we ignore the index n 
in /, so we regard / as a family of maps defined for all n. 
It is clear that our inductive definition of / on I?„ satisfies 
Eq. (fT7| . The probabilities pi,. ■ ■ ,Pn-2 G [0, 1] are free 
parameters that must be chosen to satisfy Eq. (jlSp . Note 
that these probabilities also implicitly depend on n. The 
choices of /(&) in the first two lines of the above table 
are represented by black and red nodes in the example 
shown on Fig. [2j Consider three cases: 

Case 1: a = Irt' for some t' G 23>n-2- Then f{b) = a iff 
b = llrrt' or 6 = Irt for some t S I>n-i such that f{t) ~ 
t' . These possibilities are mutually exclusive. Hence 

^ Pr[/(&) = Irt'] ^Pi+ Y. P^-[/(^) = *'] 



faeCn 



Pi 






Substituting it into Eq. (|T8)) gives a constraint 

Pi = Xn — Xn-l- 



(19) 



Case 2: a = Is'r for some s' E 'Dn-2- Then f{b) = a 
iff 6 = Is'rlr or b — Isr for some s E T>n-i such that 
/(s) = s' . These possibilities are mutually exclusive. 
Hence 



bev„ 



PT[f{b)=ls'r] = l-pn-2+ J2 Pr[/(S) = S'] 



= 1 - Pn-2 + Xn-1- 

Substituting it into Eq. (|18p gives a constraint 

Pn~2 = 1 - (Xn - -'^n-l)- 



(20) 



It says that Xn must be a non-decreasing sequence. 

Case 3: a ~ Is'rt' for some s' E Vi, t' E 'Dn-i-2, and 
i — 1, . . . ,n — 3. In other words, both s' and t' must be 
non-empty. Then f{b) = a iff 6 = Isrt' for some s E Ci+i 
such that /(s) = s', or b = Is'rt for some t E Vn-i-i such 
that fit) = t' . These possibilities are mutually exclusive. 
Hence 

^Pr[/(6) = ;s'rt']=p,+i Y. Pr[/(«) = ^'] 
bev„ seVi+i 

-i-(i-K) Y. ?'■[/(*) = *'] 

= Pi+lXi+i + (1 - Pi)Xn-i-l- 

Substituting it into Eq. ()18|) gives a constraint 



Pi+l X.i+1 + (1 - p.i)Xn-i-l = X„ 



(21) 



for each i = 1, . . . , n — 3. Let us choose 

a 4(i - 1/2) 



X- 



Ci-i 



l + l 



(22) 



Combining Eqs. (|19l20l2ip we obtain a hnear system with 
unknown variables pi, . . . , Pn-2 G [0,1]- We shah look for 
a solution {pt} having an extra symmetry 



Pi +Pn-i-i = 1 for 1 = 1,, 



,n~2. 



(23) 



One can check that the system defined by 
Eqs. (I19I20I21I23P has a solution 



Pi 



i{i + l){in - 2i - 1) 
n(n+l)(n-l) '' 



!,■ 



(24) 



Hence we have defined the desired stochastic map / : 
Vn -^ 2?n-i- This proves the induction hypothesis. 

It remains to note that for n — 1,2 the map / is 
uniquely specified by Eqs. ()17ll8p and our choice of 
Xn- Indeed, one has 2?2 = {llrr,lrlr}, Vi — {Ir}, 
and I?o = 0. To satisfy Eq. pT)) . we have to choose 
fillrr) — f{lrlr) = Ir ior n — 2 and f{lr) = for n = 1. 
It also satisfies Eq. ([T5)) since X2 — 2 and Xi = 1. This 
proves the base of induction. D 



Ground state energy for unbalanced subspaces 

Recall that the unbalanced subspace Hp.q is spanned 
by strings s € Cp^q that have p unmatched right and q 
unmatched left brackets. Our goal is to prove that the re- 
striction of H onto any subspace y-p^q with p > or g > 
has ground state energy at least n~^^^K By the symme- 
try, it suffices to consider the case p > 0. To simplify the 
analysis we shall omit the boundary term \l){l\n- Note 
that such omission can only decrease the ground state en- 
ergy. Accordingly, our simplified Hamiltonian becomes 



n-l 



(25) 



Recall that 11 is a projector onto the subspace spanned 
by states |00) - \lr), \0l) - \10), and |0r) - |rO). Let Xi{H) 
be the ground state energy of H. 

Any string s € Cp^q can be uniquely represented as 

s — UQru\rui . . . ruplv\lv2 ■ ■ ■ Ivq 

where Ui and Vj are Motzkin paths (balanced strings of 
brackets). The remaining p right and q left brackets are 
unmatched and never participate in the move 00 -f-)- Ir. It 
follows that the unmatched brackets can be regarded as 
"solid walls" that can be swapped with O's but otherwise 
do not participate in any interactions. In particular, the 
spectrum of H restricted to "Hp^q depends only on p + q 



as long as p > 0. This allows us to focus on the case 
(7 = 0, i.e. assume that all unmatched brackets are right. 
Given a string s e Cp.o, let s E {0, ^, r, a;, y}" be the 
string obtained from s by the following operations: (i) 
replace the first unmatched right bracket in s by 'x', and 
(ii) replace all other unmatched brackets in s (if any) by 
'y'. Define a new Hilbert space Tip whose basis vectors 
are |s), s € Cp^. Consider a Hamiltonian 



H=\x){xl 












(26) 



where Q^ and 6^ are projectors onto the states |0a:) — |a;0) 
and |0j/) — |yO) respectively (with a proper normaliza- 
tion). One can easily check that {s\H\t) = {s\H\t) for 
any s,t E Cp.Q. Hence the spectrum of H on Hpfl co- 
incides with the spectrum of H on Up. Furthermore, if 
we omit all the terms O^.,]^ in H, the ground state en- 
ergy can only decrease. Hence it suffices to consider a 
simplified Hamiltonian 



H^ = \x){xU + J2n,^, 



i+i 



e 



jj+i 



(27) 



which acts on Up. Note that positions of y-particles are 
integrals of motion for H^. Moreover, for fixed posi- 
tions of y-particles, any term in H^ touching a y-particle 
vanishes. Hence ff^ can be analyzed separately on each 
interval between consecutive y-particles. Since our goal 
is to get a lower bound on the ground state energy, we 
can only analyze the interval between 1 and the first y- 
particle. Equivalently, we can redefine n and focus on 
the case p = 1, q = 0, that is, assume that there is only 
one unmatched right bracket. The relevant Hilbert space 
Hi is now spanned by states 

|s) (8) |x) (Ki |i), where s e Mj-i, teMn-j- 

Recall that A^^ is the set of Motzkin paths (balanced 
strings of left and right brackets) of length k. 

We would like to treat the terms responsible for the 
motion and detection of the cc-particle as a small pertur- 
bation. To this end, choose any < e < 1 and define the 
Hamiltonian 



m 



En« 



+1 



e\x){x\ 



n-l 



i=i 



Clearly, Hf < H^, so it suffices to get a lower bound on 
the ground state energy of iff. 

Let us first find the ground subspace and the spectral 
gap of the unperturbed Hamiltonian Hq — X^^i ^jd+i- 
Note that the position of the x-particle j is an invariant 
of motion for H§ . Moreover, any projector Hi^i+i touch- 
ing the x-particle vanishes. Hence we can analyze Hq 
separately on the two disjoint intervals A = [l,j — I] and 



10 



B = [j + l,n]. It follows that the ground subspace of Hq 
is spanned by normalized states 

\^j)^\Mj-i)®\x)®\Mn-j), j = l,...,n. (28) 

The spectral gap of Hq can also be computed separately 
in A and B. Since we have already shown that the origi- 
nal Hamiltonian Eq. ^ has a polynomial gap inside the 
Motzkin subspace, we conclude that \2{Hq) > n^'-'^-^\ 

Let us now turn on the perturbation. The first- 
order effective Hamiltonian acting on the ground sub- 
space spanned by ipi, ■ ■ ■ j4'n describes a hopping of the 
x-particle on a chain of length n with a delta-like repul- 
sive potential applied at site j = 1. Parameters of the 
hopping Hamiltonian can be found by calculating the ma- 
trix elements 



ii^jl^h+ii^j) 



{^,+i\ei,+,\i^, 



Mr, 



2M„_, 



a. 



+1/ 



2Mi 



^^ 



and 



m^h+A^j^ 



1 M„_,_i M,_i 



M 



n-j 



M, 



where Mk = \M.k\ is the fc-th Motzkin number. We ar- 
rive at the effective hopping Hamiltonian acting on C", 
namely, 



Heff = |i)(i| + Erw+i, 



(29) 



where 



-aj/?,(b-)(j + i| + b- + i)01) (30) 

is a rank-1 projector. Applying the Projection Lemma 
of [isi we infer that 



Ai(7?f)>eAi(7?eff)- 



Oie'WW 



\2{H-)-2e\\V\y 



^n~l , 



where V = |a;)(a;|i -|- Yl-j=i ®i.j+i ^^ ^^^ perturbation 
operator. Since A2(i?o) — n^'^'^^K we can choose e poly- 
nomial in 1/n such that 2e||V^|| is small compared with 
^2{Hq). For this choice of e one gets 

Ai(i?f)>eAi(ifeff)-0(e')n°(^^ 

Hence it suffices to show that Ai(iJcff) > n^'~'^^\ where 
ifeff is now the single x-particle hopping Hamiltonian 
Eq. dMD- 



Let us first focus on the hopping Hamiltonian without 
the repulsive potential: 



H„ 



n-l 



This Hamiltonian is FF and its unique ground state is 

n 

Iff) - E VM,-iM^-, \j). (31) 

Our strategy will be to bound the spectral gap of Hmove 
and apply the Projection Lemma to TJcff by treating the 
repulsive potential |1)(1| as a perturbation of Hmove- 
First let us map Hmove to a stochastic matrix describ- 
ing a random walk on the interval [1, n] with the steady 
state 7r(j) = (j|.g)^- For any a,b £ [l,n] define 



P{j,k) = 6j^k - {j\Hmove\k)^ 



U{k) 



(32) 



Since •\/7r(j) is a zero eigenvector of Hmove, we infer that 
EkPij^k) = 1 and Ej<j)Pij,k) - Tr{k). A simple 
algebra shows that 

Pij^J + ^)-^^^ and P{J + 1,J)=^'-' 



2M„ 



2Mi 



are the only non-zero off-diagonal matrix elements of P. 
We shall use the following property of the Motzkin num- 
bers. 

Lemma 7. For any n> I one has 1/3 < M„/M„_|_i < 1. 
Furthermore, for large n one can use an approximation 



3" 



(33) 



where c « 1.46. 

The lemma implies that 

l<Pij,J±l)<l 

for all j. Hence the diagonal matrix elements P{j,j) are 
non-negative, that is, we indeed can regard P{j, k) as a 
transition probability from j to k. Furthermore, using 
Eq. (pij) and the above lemma we infer that the steady 
state TT is 'almost uniform', that is. 



-0(i) < 



7r(fc) 



< n°(i) for all 1 < j, fe < 



n. 



(34) 



In particular, min^ 7r(j) > 7i~ '^^ We can now easily 
bound the spectral gap of P. For example, applying the 
canonical paths theorem stated above we get 1 — A2(P) > 
l/{pl) where p is defined in Eq. ((TU)) and the canonical 
path 7(s, t) simply moves the x-particle from s to t. Since 



11 



the denominator in Eq. (|T0)) is lower bounded by n '-'(i), 
we conclude that 1 — ^2{P) > n~'~'^-^> . It shows that 

A2(i?™o.e)>n-0(l). 

To conclude the proof, it remains to apply the Pro- 
jection Lemma to H^s defined in Eq. (P^ by treating 
the repulsive potential |1)(1| as a perturbation. Now 
the effective first-order Hamiltonian will be simply a c- 
number (1|5)^ = ""(l) ^ n~^^^' which proves the bound 
Ai(ifeff)>n-o(i). 
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